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We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo 
is expected to be of a^f2 type which is likely to be relevant for the Sun, whereas the latter case refers to dynamos which are more 
likely to occur in more rapidly rotating stars. The latter exhibit oscillatory large-scale magnetic fields. We perform a parameter study 
where the shear flow is kept fixed and the rotational influence is varied in order to probe the relative importance of both modes. We find 
that in cases with shear, the streamwise component of the magnetic field is between 9 and 20 times greater than the cross-streamwise 
one, depending on the rotation rate. Oscillatory solutions are preferred either when the shear is relatively strong or weak in comparison 
to rotation, but not when the two are comparable to each other. However, exceptions to these rules also appear and in many cases 
the solution is oscillatory only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. The cases with 
rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic 
boundary conditions. The strength of total and the large-scale component of the magnetic field at the saturated state, however, are 
sensitive to the chosen boundary conditions. 
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1 Introduction 

The solar magnetic cycle is commonly thought to be a manifestation of an oscillatory large-scale dynamo 
operating within or just below the convection zone (e.g. IOssendriiveiil2003l ). A possible origin of the solar 
magnetic fields is the turbulent dynamo mechanism, where helical small-scale fluid motions and large- 
scale shear sustain the mag netic field (e.g. lMoffattl[l978l . iKrause and Radleilll980l . iRtidiger and HollerbachI 
2004 ). According to mear i-field theory, turbulen t stratified convection together with global rotatioii of th e 
Sun lead to an a effect ( Steenbeck et al. 19661 ) and large-scale differential rotation (e.g. iRiidiger'llQ^ 
Their combined effect consti tutes the aQ-dynamo which often yields oscillatory solutions (e.g. Parker 195. 
Steenbeck and Krausel[l969l ). 

However, reprodticing t he solar cycl e with direct numerical simulations still remains challenging (e.g. 
Miesch and Toomre 20091 . iKkpvla 20 111 ). Early spheri cal simu l ation s were indeed able to achieve oscillatory 
large-scale fields that propagate toward the poles ( Oilman 1983 . Glatzmaier 19851 ). Simpler Cartesian 
models with rotating stratified convection wer e less successful as only small-scale fi elds were seen (e.g . 
Nordlund et al. 1992 . Brandenburg et al. 19961 ). Only when a shear flow was added teapvla et al. 2008 



Hughes and Proctoi 



2OO2I . iKiipvia et al 



20091) tion was used (jjones and Rober ts 2000. Rotv ig and Joned 

2009bl ) . large-scale flelds were obtained. Even in the cases with imposed shear, no 
oscillatory solutions were seen although the necessary prerequisites, helical turbulence and shear 



were 



present. However, these are aspects that depend critically on the boundary conditions. Indeed. iKapyla et al. 
(l2009d ) have presented mean-field calculations of associated convection simulations that agree with each 
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other not only qualitatively in that both are non-oscillatory, but they also agree quantitatively as far as 
their excitation condition is concerned. 

Here we extend previous s tudies on large-scale dynamos due to turbulent convection in Cartesian ge- 
ometry (iKapvia et al\ 120081 . l2009bl ) to cover a larger parameter space and to explore more thoroughly 
the effects of boundary conditions on the solutions. We present runs with imposed shear and find that 
oscillatory solutions can be found in a limited part of the parameter range studied. We also report on 
rigidly rotating runs where oscillatory a^-dynamos are observed. 



2 Model 



Our model setup is the same as that of Kapyla et al. ( 20081 . 2009bl ). A rectangular portion of a star 
is modelled with a box situated at colatitude 0. The dimensions of the computational domain are 
{Lx, Ly, Lz) = (4,4, 2)d, where d is the depth of the convectively unstable layer, which is also used as 
the unit of length. The box is divided into three layers: an upper cooling layer, a convectively unstable 
layer, and a stable overshoot layer (see below). The following equations for compressible magnetohydro- 
dynamics are solved: 



VA 

'vt 



S .AyX 



{VUfA-fiovJ, 



(1) 
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^ = -^V . [/ + i V . i^VT + 2.5^ + ^ - (4) 

Vt p p p T[Z) 

where V/Vt = d/dt + {U + Uq) ■ V is the advective derivative with respect to the total (turbulent plus 
shear) flow, Uq = (0, Sx, 0) is the imposed large-scale shear flow, A is the magnetic vector potential, 
-B = V X A is the magnetic field, J = 'V x B/pQ is the current density, is the magnetic permeability, 
?7 and are the magnetic diffusivity and kinematic viscosity, respectively, K is the heat conductivity, p is 
the density, U is the velocity, g = —gz is the gravitational acceleration, and ft = r2o(— sin 0, 0, cos 0) is 
the rotation vector. The fluid obeys an ideal gas law p = pe{'y — 1), where p and e are the pressure and 
internal energy, respectively, and 7 = cp/cy = 5/3 is the ratio of specific heats at constant pressure and 
volume, respectively. The specific internal energy per unit mass is related to the temperature via e = cyT. 
The traceless rate of strain tensor S is given by 

S,^ = l{U.^^ + U,,i)-yi,V-U. (5) 

The last term of Eq. ^ describes cooling at the top of the domain. Here, t{z) is a cooling time which has 
a profile smoothly connecting the upper cooling layer and the convectively unstable layer below, where 
r-i(z) ^ 0. 

The positions of the bottom of the box, bottom and top of the convectively unstable layer, and the 
top of the box, respectively, are given by (2:1,^2,23,^4) = (— 0.85, 0, 1, 1.15)(i. Initially the stratification is 
piecewise polytropic with polytropic indices (mi, m2, ms) = (3, 1, 1), which leads to a convectively unstable 
layer above a stable layer at the bottom of the domain and an isothermal cooling layer at the top. All 
simulations with rotation use 9 = 0, corresponding to the north pole. 
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2.1 Nondimensional units and parameters 
Dimensionless quantities are obtained by setting 

d = 5 = Po = CP = /iQ = 1 , (6) 

where po is the initial density at Z2- The units of length, time, velocity, density, entropy, and magnetic 
field are 

[x]=d , [t] = ^/d/g , [U] = ^/dg , [p] = po , [s] = cp , [B] = ^ dgpofio ■ (7) 
We define the fluid and magnetic Prandtl numbers and the Rayleigh number as 

u V qd^ / 1 ds\ 

Pi. = _, Pm = -, Ra = ^ -— :r ' 8 

where xo = (PmCp) is the thermal diffusivity, and pm is the density in the middle of the convectively 
unstable layer, Zm = — Z2. The entropy gradient, measured at z^ in the non-convecting hydrostatic 
state, is given by 



cpdz J Q Hp 



where V — Vad is the superadiabatic temperature gradient with Vad = 1 — I/t, V = {dlnT/dlnp)^^, 
where Hp is the pressure scale height. The amount of stratification is determined by the parameter = 
(7 — l)eo/ gd, which is the pressure scale height at the top of the domain normalized by the depth of the 
unstable layer. We use ^0 = 0-3 in all cases, which results in a density contrast of about 23. We define the 
fluid and magnetic Reynolds numbers via 

Re = ^, Rm=^, (10) 

where /cf = 27r/d is assumed as a reasonable estimate for the wavenumber of the energy-carrying eddies. 
Note that, according to this definition, Rm is by a factor 27r smaller than the usually adopted one based 
on d instead of /cf . The amounts of shear and rotation are quantified by 

Sh=-^, Co = ^. (11) 

The denominators in Eq. (jlip give an estimate of the convective turnover time. We also use the value of 
the relative shear rate 

q = -S/no, (12) 

which is more often used in the context of disk systems, for which the local angular velocity varies like 
r^o oc r~^. The equipartition magnetic field is defined by 

Seq - {flOpUY\ (13) 



where angle brackets denote volume averaging. 
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2.2 Boundary conditions 

Stress-free boundary conditions are used for the velocity, 

U,,, = Uy^, = U,=0, (14) 

and either vertical field or perfect conductor conditions for the magnetic field, i.e. 

B^ = By = (vertical field), (15) 
Bx,z = By^z = Bz = (perfect conductor), (16) 

respectively. We may think of them as open and closed boundaries, respectively, because they either 
do or do not permit a magnetic helicity fiux. In the y and x directions we use periodic and shearing- 
periodic boundary conditions, respectively. In the runs with shear and rotation we always use vertical 
field conditions at the top and perfect conductor conditions at the bottom. The simulations have been 
made with the Pencil CODE0, which uses sixth-order explicit finite differences in space and a third order 
accurate time stepping method. 



3 Results 



We perform two sets of simulations with shear and rotation (Sets A and B) and a few exploratory runs with 
only rotation (Set C), see Tabled! In the runs with shear we keep the shear rate S constant and vary the 
rotation rate Qq. In Set A we use S = —0.05 a/ q/d and in Set B w e have S = —Q.l^J g/d. Hydrodynamical 
progenitors of the runs in Set A were used in Kapvla et al. ( 2010bl ) and the runs in Set B were obtained by 



doubling both S and VLq. In terms of q, we explore the range —10 . . . 1.99. Values of q near zero refer to runs 
with rapid (and nearly rigid) rotation, whereas small values of \q\ are associated with slow rotation. For 



and = is another special case in that a 'vorticity dynamo' (e.g. 



Elperin et al. 


2003 




Kapvla et al. 



however iKapvla et aZll2008l V We take Runs A9 and Al from lKiipvlii et a.l\ (|201ll ;) as the hydrodynamical 



progenitors for our runs in Set C with only rotation. 

In the following we discuss first the case with shear and study the behavior solutions for a range of 
values of q. We refer to these solutions as a-shear dynamos. We study separately the case without shear 
and refer to such solutions as dynamos. We focus here specifically on the case of oscillatory solutions. 



3.1 a. shear dynamos 



In our previous studies of convection-driven large-scale dynamos with shear (iKapyla et al\ l2008l . l2009bl . 
2010al V only non-oscillatory solutions were obtained. The large-scale field often had a different sign in the 
convectively unstable layer with respect to the stable overshoot layer, alth ough solutions with a single 
sign have been obtained at low i nagnetic Reynolds numbers (see Fig. 5 of Kapyla et al. 2010al ) and in 
cases where i^o = (see Fig. 7 o flKapvla et ad 120081'). Ruii s with sinusoidal she ar and rotation also show 
non-oscillatory large-scale fields ( Hughes and Proctor 2009, Kapyla et al. 2010cl ). 

As is evident form Table [H many of our current runs are also non-oscillatory; see Figure [U^a) for a 
representative result from Run A6. When q is increased from unity, the mean magnetic field begins to 
show reversals. This is particularly clear in the kinematic regime in Runs A1-A3. Cycles in the non- 
linear regime are more irregular and in some cases the dynamo mode changes to a non-oscillatory mode 
(Run Bl). With q < 0.25, another oscillatory regime is found. These runs tend to show oscillations in the 
kinematic stage but often switch to a stationary mode in the non-linear regime; see the time evolution 
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Tabic 1. Summary of the runs. In Sets A and B we use Pr — 1, Ra — 10^, Pm — 1, and grid resolution 128"^. Perfect conductor 
(vertical field) conditions for the magnetic field at the lower (upper) boundary are used. In Set C, Pm varies, Pr — 0.24, 
Ra — 2 ■ 10^, and grid resolution 256 x 128^. The boundary conditions in Set C are listed in the rightmost column of the 
table. Oscillatory (osc) and stationary (sta) modes of the dynamo are denoted in the second column from right. Question marks 
indicate that only very few sign changes are covered by the time series or irregular reversals are seen. Here, Brma — ^rma/^eq, 
where Bi-ms is the total rms magnetic field, and Bi — (S^)"'"^^/Scq. 
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Ma 


Rm 


Q 


Co 


Sh 


R 

J-'rms 


R 


R 

I3y 


Mode 


Will 111 t:llU 


Al 


0.031 


25 


1.99 


0.26 


-0.25 


1.50 


0.07 


1.38 


osc/osc? 




A2 


0.029 


23 


1.75 


0.31 


-0.27 


1.50 


0.08 


1.27 


osc/osc 




A3 


0.028 


22 


1.50 


0.38 


-0.28 


1.58 


0.10 


1.29 


osc/osc 




A4 


0.027 


22 


1.25 


0.46 


-0.29 


1.41 


0.10 


1.12 


osc?/osc? 




A4 


0.027 


22 


1.25 


0.46 


-0.29 


1.41 


0.10 


1.12 


sta/osc? 




A5 


0.028 


22 


1.00 


0.57 


-0.28 


2.65 


0.22 


2.46 


sta/sta 




A6 


0.029 


23 


0.75 


0.73 


-0.27 


3.34 


0.29 


3.16 


sta/sta 




A7 


0.030 


24 


0.50 


1.07 


-0.27 


3.92 


0.37 


3.73 


sta/sta 




A8 


0.029 


23 


0.25 


2.17 


-0.27 


4.34 


0.47 


4.08 


sta/sta 




A9 


0.026 


20 


0.10 


6.19 


-0.31 


5.70 


0.61 


5.33 


osc/sta 




AlO 


0.023 


18 


0.05 


13.7 


-0.34 


7.07 


0.75 


6.53 


osc/sta 




All 


0.011 


8 


-0.05 


-30.2 


-0.75 


- 


- 


- 


osc/— 


not run to saturation 


A12 


0.016 


13 


-0.10 


-9.87 


-0.49 


- 


- 


- 


osc/— 


not run to saturation 


A13 


0.056 


45 


-0.25 


-1.13 


-0.14 


2.21 


0.36 


1.81 


osc?/sta 




A14 


0.062 


49 


-0.50 


-0.52 


-0.13 


2.10 


0.22 


1.85 


sta/sta 




A15 


0.027 


21 


-1.00 


-0.59 


-0.30 


- 


- 


- 


sta/— 


no dynamo 


A16 


0.028 


23 


-2.50 


-0.22 


-0.28 


- 


- 


— 


osc/— 


no dynamo 


A17 


0.030 


24 


-5.00 


-0.11 


-0.26 


- 


- 


— 


osc/— 


no dynamo 


A18 


0.036 


28 


-10.0 


-0.04 


-0.22 


- 


- 


- 


osc/— 


no dynamo 


Bl 


0.090 


72 


1.99 


0.18 


-0.18 


1.15 


0.09 


1.00 


osc/sta 




B2 


0.040 


32 


1.75 


0.46 


-0.40 


2.70 


0.13 


2.54 


osc/osc? 




B3 


0.039 


31 


1.50 


0.54 


-0.40 


3.18 


0.18 


2.99 


sta/sta 




B4 


0.040 


31 


1.25 


0.64 


-0.40 


3.22 


0.20 


3.03 


sta/sta 




B5 


0.042 


33 


1.00 


0.77 


-0.38 


2.94 


0.20 


2.74 


sta/sta 




B6 


0.034 


27 


0.75 


1.25 


-0.47 


3.79 


0.34 


3.57 


sta/sta 




B7 


0.038 


31 


0.50 


1.66 


-0.42 


3.63 


0.29 


3.38 


sta/sta 




B8 


0.036 


29 


0.25 


3.50 


-0.44 


4.00 


0.33 


3.68 


osc/sta 




B9 


0.029 


23 


0.10 


11.0 


-0.55 


4.51 


0.41 


4.00 


sta/sta 




BIO 


- 


— 


0.05 


— 


- 


- 


- 


— 


— 


no convection 


Bll 


— 


- 


-0.05 


— 


- 


- 


- 


— 


— 


no convection 


B12 


0.011 


8 


-0.10 


-30.2 


-1.51 


- 


- 


- 


osc/— 


marginal dynamo 


B13 


0.018 


15 


-0.25 


-6.91 


-0.86 








osc?/— 


not run to saturation 


B14 


0.086 


69 


-0.50 


-0.74 


-0.18 


1.35 


0.19 


1.01 


osc?/sta 




B15 


0.094 


75 


-1.00 


-0.34 


-0.17 


1.16 


0.10 


0.95 


osc?/sta 




B16 


0.030 


24 


-2.50 


-0.43 


-0.54 








osc/— 


marginal dynamo 


B17 


0.032 


26 


-5.00 


-0.20 


-0.50 


1.44 


0.06 


1.20 


osc/sta 




B18 


0.035 


28 


-10.0 


-0.09 


-0.46 


1.61 


0.08 


1.35 


osc/sta 




CI 


0.021 


66 




4.60 




1.18 


0.30 


0.31 


osc / osc 


Pm = 2, pc/vf 


Clb 


0.025 


39 




3.85 




0.36 


0.09 


0.09 


osc/osc 


Pm = 1, pc/vf 


Clc 


0.032 


26 




2.96 










osc/— 


Pm = 0.5, pc/vf, dynamo decays 


Cld 


0.021 


66 




4.60 




1.18 


0.32 


0.29 


osc/osc 


Pm = 2, vf/vf 


Cle 


0.023 


72 




4.22 




0.54 


0.10 


0.09 


osc/osc 


Pm = 2, pc/pc 


C2 


0.018 


58 




17.5 




0.64 


0.10 


0.10 


osc/osc 


Pm = 2, pc/vf 



of the horizontally averaged magnetic field components from Run B8 in Figure [U^b). We find that in 
the parameter regime explored here, non-oscillatory solutions are consistently found near q = 1. This 
explains the lack of oscillatory dynamos in our previous works where we always used q = 1. This further 
i llustrates the import ance of comprehensive parameter studies instead of individual numerical experiments 
(iKapvia et ai.ll2ninbl ). 

Many runs in the q < regime, especially in Set A, are either subcritical or show very slow growth of the 
magnetic field and were not run to saturation. Convection is suppressed especially near q = due to the 
rapid rotation, decreasing the Reynolds number and thus explaining the absence of dynamo action there. 
This can seem surprising bec ause in the q < re gime the contributions to the a eff'ect due to shear and 
rotation have the same sign (iKapvla et Qi.ll2nn9bl V On the other hand, the magnetorotational instability 
can be excited for g > which may explain the more favourable dynamo excitation in that regime. 
However, we find that if a saturated dynamo is present in this regime, also the turbulence is enhanced 
(see Runs A13, A14, B14, and B15). This is associated with the generation of additional large-scale flows 
that depend on x, see Figure HI Also the large-scale magnetic fields are x-dependent. Such modes are not 
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Figure 1. Horizontally averaged horizontal components of the magnetic field from non-oscillatory Run A6 (a) and initially oscil- 
latory but ultimately stationary Run B8 (b) a-shear dynamos. 



visible in the kinematic stages of the runs. 

Even in the cases with the clearest oscillatory solutions, e.g. Run A2 in Fig. El^a), the period of the 
oscillation varies from cycle to cycle. Furthermore, the cycle period is of the order of 10^ convective 
turnover times in this run. Such a long cycle period limits the duration of the simulation to only a few 
cycles. 

The cycle frequency of a saturated a-shear dynamo under the assumption of homogeneity is given by 



eye 



(17) 



where r/x = rjt + r] is the total magnetic diffusivity, rjt is the turbulent diffusivity, and km is the 
wavenumber of the dominant mode of the magnetic field ( Blackman and Brandenburg 20021 ). Because 
this is valid in the non-liriear re gime, the quenching of r]t as a function of Rm and B can be estimated 
( Kapvla and Brandenburg 20091 ) . However, such procedure is not directly applicable here due to inhomo - 
geneity of the system and the fact that a number of other effects are also present (e.g. Kapyla et al\ 2009c ) . 
rendering Eq. (I17p inaccurate. However, provided that we are dealing with a turbulent dynamo, the cycle 
frequency is likely regulated by the value of r]t ■ This suggests that the turbulent diffusivity is quenched by 
a factor of roughly two to three in Run A2, compared to the kinematic stage of the same run. 

The phase diagram of the horizontal components of the large-scale field averaged over 0.2d < z < 0.8d 
in Run A2 are shown in Fig. [3^a). The streamwise and cross-stream field components are in antiphase in 
this case. 



3.2 Oscillatory o? dynamos 



In an earlier study w e found the appearance of large-scale magnetic fields in rigidly rotating convection 
(iKapvia ei ai.ll2009bl ). !3owever, none of the runs in that paper were run for much more than 10'^ convective 
turnover times. Although sign changes of the large-scale fields were seen (see. Fig. 7 of KKB09), the time 
series were too short to enable any firm conclusions regarding the possibly oscillatory nature of the dynamo. 

Furthermore, in similar ra pidly rotat i ng runs without inagnetic fields, the a ppearance of large-scale 
vortices has been discovered E^iT^ M. iMantere .^^ M . Here we use the hydro- 

dynamical states of runs with large-scale cyclones as initial conditions for our dynamo simulations. We find 
that a large-scale dynamo is excited provided the magnetic Reynolds number exceeds a certain critical 
value. Furthermore, as the magnetic fields become dynamically important, the cyclones decay and are 
absent in the non-linear state. The large-scale magnetic field is oscillatory in the two cases with different 
values of Co that we have considered. Figure EJ^b) shows the horizontally averaged mean magnetic fields 
from a rigidly rotating Run CI where an c? dynamo is excited. In Run CI the large-scale fields are only 
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Figure 2. Same as Fig. [l]but for oscillatory a-shcar dynamo Run A2 (a) and cP dynamo Run CI (b). 
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Figure 3. Phase diagrams for the same runs as in Fig. (2] 
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functions of z whereas in the more rapidly rotating Run C2 the large-scale fields depend also on x and y. 
Furthermore, the oscillatory nature of the solution is not so clear. Figure E^b) shows the phase diagram 
of the horizontal components of the large-scale field in Run CI. There is a phase shift of 7r/2. 

The saturation level of the dynamo is sensitive to the magnetic Reynolds number. Decreasing Rm from 
66 to 39 by doubling the value of r], decreases the saturation field strength by a factor of three (Run Clb). 
Another doubling of rj shuts the dynamo off (Run Clc). 

Our standard setup in the present paper is to use perfect conductor boundaries at the bottom and 
vertical field conditions at the top. Changing also the lower boundary to vertical field conditions produces 
no discernible difference in the solution (Run Cld). However, imposing perfect conductor conditions on 
both boundaries decreases the saturation strength to less than a half of the standard setup and decreases 
the fraction of the large-scale field (Run Cle). We have not, however, studied the Rm-dependence of the 
saturation field strength in this case (cf. iKapyla et al. 2010al ) . 



4 Conclusions 



We have presented results from simulations of turbulent magnetized convection both with an imposed 
shear flow using the shearing box approximation (Sets A and B) and in rigidly rotating cases (Set C). In 
accordance with previous results, we find the generation of dynamically important large-scale magnetic 
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Figure 4. Mean velocities (upper row) and magnetic fields (lower row) averaged over the y-direction and time from Run B14. The 
velocities are shown in units of s/dg and the magnetic fields in units of volume averaged equipartition field Beq ■ 



fields. In distinction to our earlier studies, we here vary the relative importance of rotation with respect 
to shear by covering the range q = —10. . . 1.99 of the relative shear rate q. We find that for q = 1 the 
solutions are always stationary, which is in accordance with earlier results. For large q, i.e. slower rotation, 
and for g <C 1, oscillatory solutions are found. These trends are particularly clear in the kinematic regime. 
In the saturated state, however, we often find that the dynamo switches from oscillatory to stationary. 

In the rigidly rotating cases of Set C, all the dynamo solutions are found to be oscillatory. Large-scale 
vortices, present in the hydrodynamic state, are no longer found in the saturated state of the dynamo. 
Usage of a perfect conductor boundary condition instead of a vertical field condition, allowing for magnetic 
helicity fluxes, is found to decrease both the total saturation field strength and the strength of the large- 
scale field with respect to the total magnetic field. 
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